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(AUV) is addressed. The highly nonlinear nature of the vehicle 
dynamics and the requirement for fast maneuvering call for 
robust control techniques. In particular, Variable Structure 
Control ( VSC ) combined with Adaptive Control (AC) techniques 
seem to yield satisfactory performance in terms of robustness, 
capability to adjust to different operating conditions, and 
speed of response. Also, linear robust techniques based on LQG 
and robust observers are presented to address the case when 
the whole state (in terms of pitch rate, pitch and depth) is 
not available for measurement. 
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I. INTRODUCTION 



In the last few years, considerable interest in Autonomous 
Underwater Vehicles (AUV) has arisen. Possible utilization 
scenarios for an untethered submersible have also grown and 
include: sea floor mapping, target identification, remote 
reconnaissance, object recovery, etc. A typical mission would 
include downloading of the mission program, remote positioning 
and launch, subsequent vehicle recovery, and data collection. 
AUV vehicles now in design or production can employ basically 
any shape from a box shape suited to offshore work to a body 
of revolution suited for high speed maneuverability or long 
distance missions. The AUV group at the Naval Postgraduate 
School (NPS ) has chosen the shape shown in Figure 1 which was 
derived from the Swimmer Delivery Vehicle (SDV) [Ref. 1]. Part 
of the rationale was that, since extensive data exists on that 
particular design, it would facilitate calculation of a 
reasonably accurate analytical model for a similar but smaller 
scale model. This is particularly advantageous to the design 
of a control system for the AUV, since knowledge of an 
accurate dynamic model is in general the basis of a reliable 
control design. Also, tests conducted by computer 
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simulations yield results which are closer to the physical 
performance of the vehicle. 

The initial goal of the AUV group is to build a working 
prototype AUV that will be able to expand its capabilities to 
meet the increased complexity of future missions. 




As with all analytical models of submersibles , the 
equations of motion are highly nonlinear and vary with speed 
and the amount of equipment onboard. Adding to that is the 
difficulty in obtaining exact dynamic coefficients. Because 
of the high level of nonlinearities and model uncertainties, 
we have to consider more sophisticated design techniques. 
Recently, an adaptive controller has been designed by Schwartz 
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[Ref. 2] where the parameters of a compensator are adjusted 
on-line on the basis of an estimated transfer function. This 
approach has not only provided insight for the overall control 
design, but it has also provided parameters of the linearized 
dynamics . 

The adaptive controller by Schwartz has proved to yield 
satisfactory behavior in terms of speed of response and 
robustness in the presence of model changes. On the other 
hand, the computational complexity involved in the adaptive 
controller calls for investigation of simpler techniques which 
do not require complex manipulation of matrices to be 
performed in real time. 

With these motivations in mind (mainly adaptability and 
robustness in the presence of uncertainties) , we have 
investigated an alternative control technique based on 
Variable Structure Control (VSC) . It will be shown that a 
simple and robust nonlinear controller can be designed, 
combined with an adaptive loop which compensates for some of 
the model uncertainties. 

The exciting part of the VSC approach is that the design 
can be based on a nominal model (usually available from test 
runs and/or physical insight) and a bound on the uncertainty 
of the model. This gives to the designer the possibility of 
using this simple nominal model, possibly linear, while still 
preserving the global stability of the controlled system. 
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Situations which exhibit coupled dynamics can be handled by 
VSC design techniques, and simple and robust controllers can 
be designed by including the effects of coupling in the model 
uncertainty. 

A further aspect of this research is the investigation of 
the performance of a robust linear controller in the presence 
of uncertainties and nonlinearities. In particular we address 
the problem of control when the full state is not available 
for measurement. This is the case, for example, when in a 
diving maneuver we measure depth alone, rather then depth, 
pitch and pitch rate. It turns out that the design of a 
"robust" observer (in a sense which will be clarified in a 
later section) together with an optimal controller can still 
be satisfactory, under limited ranges of the operating 
conditions. 

An outline of the thesis is as follows. Chapter II covers 
modeling of the AUV and stability analysis, while Chapter III 
addresses the controller design by Variable Structure Control. 
In Chapter IV a controller design based on the theory of 
Robust Observers is presented with applications to the AUV. 
Alternative hardware realizations are surveyed in Chapter V, 
while conclusions and recommendations for future research are 
the subjects of the concluding chapter. 
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II. MODELING OF THE AUV 



A. BASIC DEFINITIONS 

The dynamics of underwater vehicles are nonlinear and are 
affected by uncertainties due to modeling errors (data 
sampling, sensors, electrical noise, physical measurements, 
hydrodynamic coefficients, etc.) and external disturbances 
(currents, tether, surfaces, etc.). In general, a linear 
controller for a vehicle which moves in "n" degrees of freedom 
should be linearized about its respective axes at sufficiently 
different speeds and then some form of "gain scheduling" is 
usually implemented. A mathematical model of the AUV can be 
obtained by accounting for the hydrodynamic forces acting upon 
it. In particular, the model results in a set of highly 
nonlinear differential equations in state space form 

x = f (x,u, t) (2 . la) 

where x is the vector of position and velocities (both 
absolute and angular) of the vehicle. The state vector x is 
comprised of 

x = surge <f> = roll 

y = sway B = pitch 

z = heave = yaw 

which represent position (absolute and angular) in a global 
coordinate frame, and 
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u = surge rate p = roll rate 

v = sway rate q = pitch rate 

w = heave rate r = yaw rate 

which represent velocities in a body fixed coordinate frame. 
Therefore the full state vector x in the differential equation 
becomes x=[u,v,w,p,q,r,x,y,z,^,0 ,V>] 1 • 

In a diving maneuver the only states of interest for this 
diving controller are depth (z) , pitch (0) and pitch rate (q) 
while the input we consider is the dive fin (stern plane) 
deflection angle (d) . Therefore the dynamic equations become 
x=f (x, d) (2.1b) 

where 

x=[q,0,z] (2.1c) 

and d is the stern plane dive command. Simulations conducted 
using the full nonlinear model [Ref. 2] show that, within the 
neighborhood of operating conditions (i.e., cruising at a 
constant speed) , a linear model yields a fairly reasonable 
fit. This has also been confirmed by data collected from the 
prototype [Ref. 1]. 

The basic structure of the dynamic model considered is 
shown in Figure 2. The nonlinear model accounts for 
perturbations around the linear model which are caused by 
various hydrodynamic effects (currents, cavitation, surface 
effects, etc.). These effects have the tendency to perturb 
the AUV from the nominal trajectory, thus acting like an extra 
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Figure 2 AUV Dynamic Model 

component on the stern plane. Therefore, we can model this 
effect by the feedback signal fb 2 (t) (shown in Figure 2) which 
adds a perturbation to the dive plane as 

fb 2 (t) = f (q(t) ,©(t) ,z(t) ,t) . (2.2) 
We can assume that the perturbation is time-varying and 
depends on pitch rate, pitch, and heave (depth) . 

B. LINEAR MODEL 

To develop a specific linear model which approximates the 
dynamics, let us consider a structure of the form shown in 
Figure 3(a). For simplicity, let G(s) be a first order 
transfer function and also note that this model does not 
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contain a feedback from ©(t) to the input (i.e., the gain K 0 
in Figure 2) . There are several reasons to believe that this 
feedback gain is required. First, physically, the equalizing 
moments are 

Jq f F,(dCt)) + F 2 (©(t)) + F 3 (2.3) 
where F.,(d(t)) is the torque around the center of gravity due 
to the stern plane angle, F 2 (©(t)) is the torque due to the 
displacement between the Center of Gravity (CG) and the Center 
of Buoyancy (CB) , and F 3 is the balance of the hydrodynamic 
moments acting over the body of the vehicle (including 
hydrodynamic added inertia, drag, and angle of attack induced 
moments) . For small values of 0(t) , this effect is 
proportional to 0(t) itself, as diagrammed in Figure 3(b). 
Second, it has been shown in Reference 2 that the response of 
the dynamics relating d(t) and 0(t) from an imparted impulse 
to the stern plane yields a pitch 0(t) which decays to zero. 
This indicates that all the poles of the transfer function 
between d(t) and 0(t) are stable. Finally, as speed increases, 
we expect that the stabilizing effect due to the feedback gain 
K 0 decreases because other hydrodynamic forces become more 
dominant. The result is that the dominant poles of the 
transfer function from d(t) to ©(t) follow a root locus path 
as shown in Figure 3(c). [Ref. 2] 

Following this formulation, the linearized state-space 
mathematical model for the AUV in a dive maneuver is 
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- Vz = V sin 0 

(using small angle 
approximations) 

VZ = -V 0 
Z(s) = -(V/s) 0 (s) 





Figure 3 AUV Characteristics 
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(2.4a) 



z = [ 0 0 -V ] x (2.4b) 

x = [ q 0 z ] » (2.4c) 

where we define the "normalized depth" as 

z = ( — 1/v ) z (2 . 4d) 

with v being the forward velocity of the AUV. The fact that 
the system is in controllable canonical form will be 
instrumental in the design of the control system. The 
parameters a, c and K<, in (2.4a) depend on the operating 
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conditions (speed, primarily) of the vehicle, and can be 
determined on the basis of test runs of the AUV. On the basis 
of these runs we can determine nominal values a 0 , c 0 and K 0 for 
the parameters and write the dynamics as 

x = A 0 x + b (d +f) (2.5) 
and design the controller based on this model. In the 
following chapters control system design techniques based on 
this model will be presented. 

C. STABILITY METHODOLOGY 

Most methods of stability criteria, including Routh and 
Nyquist, are not applicable to nonlinear systems. The second 
method of Lyapunov is the most general test for stability 
analysis. According to this approach we define stability (in 
the sense of Lyapunov) and a criterion for stability as 
follows . 

Definition: Given a nonlinear system 

x = f (x,u, t) (2.6) 
and an equilibrium point x e , such that f(x e ,u,t) = 0, we say 
that x e is a stable equilibrium point (in the sense of 
Lyapunov) provided that for all e > 0 there exists a S £ > 0 
(where 6 depends on e) such that all trajectories of (2.6) for 
which ||x(0)-x e || < e do not leave a ball of radius 6 e around x e , 
i.e., | x ( t ) — x e || < 8 £ for all t > 0. 
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Definition: If x e is stable and 

lira x(t) = x e (2.7) 
for any initial condition x(0) , then x e is an asymptotically 
stable equilibrium point. [Ref. 3] 

A method to determine whether an equilibrium point is 
stable or not is by the second method of Lyapunov. This method 
is quite convenient for the stability analysis of nonlinear 
systems since it does not require explicit solutions of the 
differential equations. To briefly summarize the method, 
consider a system described by the differential equation 
x=f(x,t) with f ( 0 , t ) =0 for all time. If there exists a 
differentiable scalar function V(x,t) such that V(x,t) is 
positive definite and V(x,t) is negative definite, then the 
equilibrium point is asymptotically stable. A precise 
statement of the method is presented in Reference 3 . 

Given this definition of stability and the AUV model 
presented earlier, we were ready to design a controller. A 
nonlinear design method using Variable Structure Control (VSC) 
was then followed. 
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III. VARIABLE STRUCTURE CONTROL 



A. INTRODUCTION 

The nonlinear dynamics of the AUV can be described by 

x = f(x,u,t) (3.1) 
and present a high degree of uncertainty due to several 
factors. In spite of the highly nonlinear behavior we can 
still approximate the dynamics of the AUV by a linear 
component and a nonlinear perturbation 

x = Ax + b(d + f (x) ) . (3.2) 
In the particular case of a diving maneuver, the state of 
equation (3.2) is given by x = [q,0,z"] ' and the dynamics (3.2) 
assume a controllable canonical form structure (see (2.4a)). 

The aim of the controller in a diving maneuver (like in 
many similar problems) is to drive the state vector x to track 
a desired state x d . In other words, given a desired state Xd 
G of the form 

Xd(t) = [x d °(t),... / X d (n ’ 1) (t)]' (3.3) 
we want to design a controller such that the error signal 

e(t) = x ( t ) - Xd(t) (3.4) 
tends to zero as t -* oo. Furthermore we want all signals in 
the loop to be bounded for all time. 

In our particular case the desired signal x d is 
represented by the desired depth, pitch and pitch rate. Since 
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in steady state (i.e., cruising at a constant depth) pitch and 
pitch rate are zero, we define x d as 

x d = [0 0 -Zd/v] ' (3.5) 

with z d being desired depth. 

Due to the presence of the uncertainty f(x) in (3.2), a 
control solution which guarantees a sufficient stabilty margin 
has to be devised. The Variable Structure Control (VSC) (also 
called the Sliding Mode (SM) ) technique is at the basis of the 
controller we propose for the AUV. 

B. VSC THEORY 

The idea behind using the VSC method on a dynamic system 
x = f (x,u, t) (3.6) 

with x e vC and u 6 S 1 is to drive the state from any initial 
condition x(0) onto a sliding surface 

s(x) = 0 (3.7) 

of the state space in a finite time and to keep the state on 

the surface for all subsequent times. Although the surface 
(3.7) can be chosen in a fairly arbitrary fashion, it 
nevertheless must be associated with stable dynamics. In 

particular, the surface s(x(t)) = 0 for all t > t 0 must imply 
lim x(t) = 0 (3.8) 

In most applications the sliding surface is a linear function 
of the state, in the sense that (3.7) is chosen as 

s (x) = c'x (3.9) 

for some vector c' e vC. 
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From the statements presented in the last paragraph, the 
terminology of Sliding Mode is evident, since the state, once 
taken from the initial condition x(0) onto the surface s(x) 
= 0, "slides" to zero on this surface by virtue of (3.8) (see 
Figure 4) . Apart from the requirement of being associated with 
stable dynamics, a surface in $0 is a sliding surface provided 



X(0) 


k 

\/ 




v N A 

^^^-S(X)=0 



Figure 4 Sliding Surface 



that we can determine a control input signal such that 

s { X ( t ) }s{x(t) } < 0 (3.10) 

whenever s{x(t) ) ^ 0. The reason behind condition (3.10) is 
that the definition of the Lyapunov function 

V (x) = 0 . 5s 2 (x) (3.11) 



yields 



V(x) 



s(x)s(x) 



(3.12) 
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along the trajectories of (3.6), and therefore the condition 
V (x) < 0 (given by (3.10)) makes s 2 (x) a monotonically 
decreasing function so that s{x(t) } tends to zero. Moreover, 
if we "strengthen" condition (3.10) by imposing 

s{x(t) }s{x(t) } < -n |s{x(t)}| (3.13) 
with rj a positive constant, then it is easy to see that the 
sliding surface s (x) = 0 is reached in a finite time. This 
can be seen by writing (3.13) as 

s(x(t)} < -T?sgn{s (x(t) } (3.14) 
where we define the sgn function as 



sgn (s) 



+ 1 
-1 
0 



if s > 0 
if s < 0 
if s = 0 



(3.15) 



On the basis of these conditions, if we restrict ourselves 

to linear sliding surfaces (s(x) = c'x as in (3.9)), then we 

obtain from (3.14) and (3.6) 

s{x(t) } = c'f(x,u,t) < -r?sgn{s(x)} (3.16) 

for which the control u(t) must be such that 

u + (t) if sgn{s(x)} > 0 
u (t) = (3.17) 



u.(t) if sgn{s(x)} < 0 
and u + and u. are such that 

c 1 f ( x , u+ , t ) < -T) (3.18a) 

C'f(x,u.,t) > r? (3.18b) 

with r? being a positive constant parameter of the controller. 
As a matter of fact, rj does not have to be a constant and, in 
general, it can be made a function of the state r? (x) • From the 
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definition of the control action (3.17) the terminology of 
Variable Structure Control (VSC) becomes evident, since the 
control assumes different structures (u + and u.) according to 
which side of the sliding surface the state is on (i.e., 
according to the sign of s(x)). 

C. APPLICATION OF VSC TO TRACKING 

A particularly interesting situation arises when the 
dynamic model (3.6) can be written in the form 

x = f (x) + b (x) u (t) (3.19) 
where nominal values of f and b corresponding to f and b in 
(3.19) are known. For simplicity assume that b(x) is 
completely known (this will be relaxed at the end of the 
section). We can write f(x) as 

f (x) = f(x) + Af (x) (3.20) 
with Af representing the uncertainty in f(x), and we assume 
to know some bound on |Af(x) | to be specified later. Given a 
desired state x d (t) , we want to apply the sliding mode 
technique in order to move the state x(t-) of the system to 
track the desired trajectory x d (t) . In order to do this we 
define the error as 

e(t) = x ( t ) - x d (t) (3.21) 
and a generic surface as 

s (e) = c'e (3.22) 
for some vector c' e ’SC . Combining (3.19), (3.20), (3.21) and 
(3.22) we obtain 
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c' (x - x d ) = c'f(x) + c'Af(x) 

-c'x d (t) +c'b(x)u(t). (3.23) 

An important feature of (3.23) is that the right hand side is 
the sum of terms which are known at all times (c'f(x) 
c'x^t)}, an uncertain term (c'Af(x)} and the control term 
(c 'b (x) u (t) } . For this reason we can separate the control 
input u(t) into two terms as 

u ( t ) = u(t) + u(t) (3.24) 

with u(t) compensating for the "nominal" dynamics, i.e., 

c'f(x) - c'x d (t) + c'b(x)u(t) ~ 0 (3.25) 

and u(t) compensating for the uncertainty Af. Combining 
(3.23), (3.24) and (3.25), we obtain 

c ' e = c ' Af + c ' bu (3.26) 

and we can determine the input component u in order to drive 
the state e(t) on to the sliding surface as 

s { e (t) } = c'e(t) = 0. (3.27) 

This determination can be accomplished by choosing u(t) so as 
to satisfy (3.18), which leads to the control 

u(t) = -F(x,t) sgn{c'e(t) } (3.28) 

with F(x,t) a known function of the state, such that 

F (x, t) > |c' Af (x) |/| c'b(x) | . (3.29) 

Clearly a condition which has to be satisfied is that 
c'b(x) +■ 0 for all time. 
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In conclusion, if c' is such that c'e(t) = 0 (which 

implies e(t) -► 0) , then we can say that the control input from 
(3.28), (3.29) and (3.25) 

u ( t ) = u ( t ) + u(t) = -(l/c'b) {c' f (x) - c'x d (t)} 

+ F (x, t) sgn (c 1 e (t) ) (3.30) 

yields the following: 

• the error e(t) is driven to the sliding surface 
c'e = 0 in a finite time, and 

• e (t) — * 0 as t — » oo. 

A very important point is that, once the state is on the 
sliding surface, the decay of e(t) is determined by the vector 
c' only . [Ref. 4] 

D. VSC EXAMPLE 

To illustrate the sliding mode design process, let us 
implement it on a second order system described by 

x (2) = f(x,t) + b (x) u (t) + p(t) (3.31) 

where we define the state x(t) = [x (1) x] ' , the desired state 
x d = [ x d 0) x dl ' and the vector c' = [1 A] (with A > 0) . While 
f(x,t) is not known precisely, we assume that nominal values 
f(x,t) are known, for which 

f(x,t) = f(x,t) + Af ( x , t ) (3.32) 

with Af ( x , t ) having a known upper bound of 

F (x, t) > | Af (x,t) | . (3.33) 

The estimate f(x,t) is available from several sources where, 
for example, a nominal model could be generated from 
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experimental data and parameter estimation techniques. If the 
control gain b(x,t) is also uncertain, then we assume that it 
is bounded by a function B(x,t) where the gain is known within 
a certain ratio 

1/B (x, t) < b(x,t)/b(x,t) < B (x, t) . (3.34) 

Also, the perturbations p(t) are assumed to be unknown and 
bounded by a continuous time function 

P(t) > |p(t) | . (3.35) 

The dynamics of s(x,t) are required so (3.9) can be 
differentiated with respect to time to get 

s(x,t) = x - x d + Ae. (3.36a) 

Then by substituting (3.31) into (3.36a) we obtain 

s(x,t) = f(x,t) + b(x,t)u(,t) + p(t) - x d + Ae (3.36b) 
where f(x,t) and b(x,t) are as described in (3.32) and (3.34) 
respectively. The system is stable if s(x,t) converges to 
zero, which can be assured if u(x,t) is chosen to satisfy 
(3.13). The total control consists of two parts (as in 
(3.24)), one for the known or estimated part of the dynamics 
(u) and the other for the uncertain or nonlinear part (u) . 
Assuming initially that there is no uncertainty in either the 
control matrix (b(x,t) = b(x,t)) or the transition matrix 

(f(x,t) = f(x,t)) and that there is no external perturbation 
present (p(t) = 0) , then, for this example, the known control 
is obtained by setting (3.36b) to zero so that 

u(x,t) = -(l/b(x,t) ) [f (x,t) + \e - x d ] . (3.37) 
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This control law does not yet satisfy the stability criterion 
because f (x , t ) fi (x, t) and the perturbations p(t) have not been 
accounted for. The complete control law will have another term 
(u) that is discontinuous across the surface in order to 
satisfy ss<-r?|s|. Now with all the uncertainties included, 
where s=s(X,t), s becomes 

s = f+ Af + bu + p-x d +Ae (3.38) 

and, after substituting u=u+u into ss, it becomes 

ss = [Af + p - u]s < -r?|s|. (3.39) 

This results in u=-k(x,t) sgn(s) so that 

u(x,t) = u(x,t) - k(x,t)sgn(s) (3.40) 

where k(x,t) is determined from the bounds on the 
uncertainties and perturbations previously estimated as 

k (x, t) = [ F (x , t ) + P(t) + r/] . (3.41) 

This variable k(x,t) is the mechanism by which the system 
uncertainties are accounted for and which will cause the 
discontinuous part of the control to compensate for their 
effect. This also insures that s 2 (x,t) is a Lyapunov function 
which in turn guarantees stability. 

If the control gain is also uncertain, then the control 
must be changed to 

u(x,t) = [u(x,t) - k(x,t) sgn(s) ]/b(x, t) (3.42) 

where the discontinuous term now is 

k (x, t) = B [ F (x , t ) + P(t)] + (B-l) | U (x , t) | (3.43) 
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and B is as defined in (3.34) [Ref. 5]. Again notice that the 
control discontinuity has become larger to overcome the 
increased uncertainty in the system. 

This discontinuity causes a chattering problem which is 
present due to the abrupt nature of the sign function, so a 
smoothing operation needs to be implemented. This can be done 
by replacing the sign function with a saturating function 
which has a user defined slope or boundary layer between the 
upper and lower limits of plus and minus one. This provides 
a linear region around the surface boundary and smooths the 
control action (see Figure 5a and 5b) . Now the designer has 
another variable to use, <f> (called the boundary layer 
thickness and whose magnitude depends on the level of system 
uncertainties) , along with X to incorporate into his design. 
These variables can be made time-varying to account for times 
when the uncertainties increase or decrease, thus increasing 
overall performance. [Ref. 6] 

The questions of rigorous mathematical uniqueness, 
existence, and continuability have been addressed in Reference 
5 and Reference 7. Though bang-bang controllers are 
attractively simple and optimal control specialists have 
pretty much perfected their use, the problem is that the 
differential equations governing their use have 
discontinuities on the right hand sides. Therefore, normal 
uniqueness and existence theory for differential equations no 
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longer applies. Filippov proves the existence and 
continuability for his own solution concept [Ref. 7]. For the 
question of uniqueness, it can be shown that, by satisfying 
(3.13), only one solution to the discontinuous differential 
equation can apply at any given time [Ref. 5]. This is 
equivalent to saying that the derivative of the state 
trajectory must always point towards the sliding surface. 
Although VSC in itself is a powerful design methodology, 
incorporating an adaptive portion can increase the robustness 
and performance of the controller. 
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E. VSC WITH ADAPTIVE COMPENSATION 



With the linear model developed in Chapter II, we will 
now determine a simple adaptive controller capable of tracking 
the state of a reference model in the presence of 
uncertainties in the AUV dynamic parameters. The state space 
description of this model is 

x ( t ) = Ax (t) + bd (t) + f(x) (3.44) 
where d(t) is the control input and f(x) is the nonlinearity 
which comprises model inaccuracies, external perturbations 
and changing parameters. This depth controller design is based 
on the model dynamics with an adaptive loop to place the 
eigenvalues for the desired linear part and a switching input 
based on VSC design for the nonlinear part f(x). The 
controller should drive the AUV to a desired depth z d and keep 
it there, despite system uncertainties and external 
disturbances. The difference between the actual state 
[q, ©, z/ (-v) ] ' and the desired state [ 0 , 0 , z d / (-v) ] 1 (where v 
is the vehicle speed) is defined as the error state 

e(t) = [q,0, (z-z d )/(-v)] '. (3.45) 
Due to the presence of integral action in the vehicle model 
(3.44), we can include the constant desired depth in the 
initial condition of the integrators and write (3.44) with e 
as the state (rather than x) • 

Now let us choose a model matrix A m with eigenvalues 
determined by the desired closed loop response. Clearly Am 
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has eigenvalues in the stable region. If we choose A m such 
that the pair (A^b) is in controllable canonical form, then 
it is a simple exercise to show that the vehicle dynamics 
(3.44) can be written as 

e (t) = A^Ct) + b (d (t) + K'e(t)) + f(e) (3.46) 
where the vehicle dynamics (at the current operating 
conditions) determine the required gain K. This dynamic model 
(3.46) is at the basis of an adaptive controller which will 
drive the error state e(t) to zero and track the ordered 
depth. The surface s(e) = 0 is defined (as described in (3.9)) 
as a linear combination of the error states, namely 

s(e)=c'e(t). (3.47a) 
For this particular formulation, let c' be the left 
eigenvector of the matrix A,,, associated to any of its stable 
eigenvalues, -A(i.e., c'A,,, = -Ac'). By the fact that A m is in 
companion form by assumption, it is possible to show [Ref. 8] 
that the entries of the vector c* are the coefficients of the 
polynomial having all other eigenvalues of A m as roots. Since 
q=0, 9=z and z d is constant, we can write s(e) as 

s(e) = c 2 e (2) + c^ 15 + c 0 e (3.47b) 
where e= (z-z d )/ (-v) , and c'=[c 2 c, c 0 ] . All this implies that 
the surface s(e)=0 satisfies one of the requirements of being 
a sliding surface, i.e., 

s (e) = 0 =* lim e (t) = 0 (3.47c) 
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On the basis of this definition we can determine an adaptive 
controller which drives the state error e onto the surface 
s(e)=0 as follows. 

In the model (3.46) let us assume that bounds 

K, min < K,- < K, max (3.48) 

on the entries of the vector K, and a bound 

F(e) > | c ' f (e) /c *b | (3.49) 

on the nonlinearities are known to the designer. Then, under 
this assumption, the control input is 

d(t) * -K' (t)e(t) + F (e) sgn{ s (t) } (3.50) 

with K the adaptive gains defined as 

K(t) = -ct {K ( t) } - M e(t)s(t). (3.51) 

The expression <7 r {K(t) } is given by 

0 if Kj min < Ki < K, MAX 

a f ( K ( t ) ) = -a(K i (t) - K^ m ’ n ) if K f (t) < K, m ’ n 

-a(K, (t) - K, max ) if Kj (t) > K,- MAX (3.52) 

(where a is a positive constant), hence (3.51) yields a closed 
loop response which is exponentially stable and e(t) — ► 0. 
[Ref. 4] 

Proof : Using the fact that s=c * e and c'A m =-Ac', we can 
write (3.46) from (3.47) and (3.50) as 

s(t) + As(t) = c'bK'(t)e(t) + c'f(e) 

-c'bF(e) sgn{s(t) } , (3.53) 

where K = K - K is the parameter error. Define the Lyapunov 
function as 

V(s,K) = 0 . 5 ( S 2 + gK'K) (3.54) 
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where g = c'b/n (a positive quantity). The time derivative of 
(3.54) is 

V(s,K) = -As 2 (t) - c'bF(e) | s (t) | + c'f(e)s(t)) 

-gK’(t)a{K(t)}. (3.55) 

Now we can see that the term K' (t)a(k(t) } is always 
nonnegative by definition of the function £(K(t) = 
[(7 1 (k(t) ) , • • • f (7 n (K(t) ) ] ' in (3.52). Also, by the definition 
of the bound F(e) on the nonlinearity, the term (c'bF(e) |s(t) | 
- c'f(e)s(t)} is always nonnegative. Therefore 

‘ V(s,K) < -As 2 (e (t) ) < 0 (3.56) 

along the trajectories of the system. This implies that 
s(e(t)) is always bounded and also that the adaptive gains 
K(t) are bounded. This, combined with (3.56), yields 

>J + Vdt < -J + ~ V(t) dt < V(0)-V(oo) < oo (3.57) 
from which we deduce that 

lira s (e (t) ) = 0 (3.58) 

and e(t) -> 0, which proves the result. 

F. SIMULATIONS AND RESULTS 

The controller (3.50) has been implemented in MATLAB and 
the basic flow chart is outlined in Figure 6. The performance 
is satisfactory at all speeds and, depending on the desired 
closed loop eigenvalues (i.e., the choice of the matrix A„,) , 
the response could be tailored towards minimum rise time or 
minimal overshoot. While we used a saturation function instead 
of a sign function in our program, our process achieves 
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VSC Algorithm 




Figure 6 VSC Algorithm Flow Chart 

further smoothing of the discontinuous boundary by use of the 
integral process in determining the feedback gains Kj . The 
simulation results are shown in Figure 7 while Figure 8 shows 
a plot of typical feedback gains. 

These results show that the controller performs 
satisfactorily over a wide range of operating conditions, 
ranging from 100 to 500 rpm of the thrusters. The different 
rise times appearing in the depth plot in Figure 7 are 
consistent with the different speeds of the vehicle. In all 
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Figure 7 VSC Performance Plots 

three runs the dive plane action (bounded within ± 0.4 
radians, or ± 23 degrees due to the physical constraints) does 
not show any chattering due to the linear saturation adopted 
for the controller (3.50), rather than the sgn function 
formulated in the theory. The result is a smooth operation of 
the dive fin (stern plane) . 

The adaptive gains K-| , K 2 and K 3 plotted in Figure 8 stay 
within reasonable values and converge to constants which 
depend on the operating conditions. The surface S is also 
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Figure 8 VSC Gains 

plotted and is shown to be driven to zero and kept there. For 
these particular plots the speed of 300 rpm is assumed. 

Considering the simplicity of the controller, we could 
conclude that this control scheme is satifactory for on-line 
implementation on the AUV under development at NPS . Although 
the VSC method is an efficient nonlinear design method, other 
linear design schemes are available which can also satisfy the 
robustness requirement. 
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IV. LINEAR ROBUST CONTROL OF THE AUV 



A. INTRODUCTION 

The VSC described in the previous section has been shown 
to exhibit satisfactory robustness properties in the presence 
of nonlinearities and unmodeled dynamics. A major requirement, 
however, of VSC is that the states must be available for 
measurement. For the AUV on a dive maneuver this is not much 
of a drawback since the state signals (pitch rate, pitch and, 
depth) are provided by gyros and the depth cell. 

There are situations, however, in which we might want to 
be able to control the vehicle with incomplete state 
information. This is the case when it is important to provide 
for reliability in the sense of ensuring the success of the 
mission, even in the presence of failure of a gyro or its 
circuitry. In these cases a model based on an observer 
provides for the missing measurements by estimating them from 
the other available signals. 

It is a well known fact that robust design techniques 
based on full state feedback are bound to lose their 
robustness properities when the state is replaced by an 
estimate from an observer. This has been pointed out by Doyle- 
Stein [Ref. 9] in the context of robust linear controllers. 
In the next section we introduce the notion of robust 
observers designed to preserve the robustness properties of 
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the compensator. The performance of the AUV with this type of 
observer is given in the last section. 

B. ROBUST OBSERVERS AND THE DOYLE-STEIN CONDITION 

The design of a compensator usually starts with a system 
description such as 

x = Ax + Bu (4.1a) 

y = C'x + Du (4.1b) 

where A is the plant transition matrix, B is the control 
matrix, C is the observation matrix, and D is the feedthrough 
matrix (usually it is zero) . Because of the Separation 
Principle, the controller and observer can be designed 
independently. Designing the controller is usually the first 
step and starts with finding the gains of the feedback control 
law (u=-Gx) where the required response is determined from 
system specifications. The gains are then selected by a pole- 
placement formula such as the Bass-Gura method or Linear 
Quadratic Regulator (LQR) methods, based on a quadratic 
performance index, and hence solve a Riccati equation 
associated with the system (4.1) . Real life constraints on the 
physical system limitations (such as the power supply) limit 
the size of the control u, and this limitation can easily be 
imbedded in an LQR design. 

When all state signals are not available for measurement, 
an observer provides for their estimates from input/output 
measurements. Using a Luenberger observer of the form 
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X = AX - Bu + K(y-c'x) (4.2) 
it is well known that x-x -» 0, provided A-Kc' has stable 
eigenvalues. With the observer, the control input becomes 
u=-Gx. The use of an estimated state (rather then the state 
itself) might seriously affect the robustness of the 
controller. Even the seemingly sensible solution of making 
the estimated state x converge to the actual state x as fast 
as possible (i.e., decrease the negative real part of the 
eigenvalues of A-Kc') does not, in general, improve robustness 
of the closed loop system [Ref. 9] . 

A significant problem concerning a lack of robustness 
(i.e., sufficiency of the gain and phase margins in the 
presence of parameter variation) has been associated with 
systems whose observers were designed solely on the assumed 
sensor noise parameters [Ref. 10] . Therefore, the issue of 
robustness of the closed loop dynamics should be investigated 
as part of the design process of an observer. A typical LQR 
controller using full state feedback can have gain and phase 
margins well in excess of six decibels and 60 degrees 
respectively; however, using an observer, the margins can be 
reduced significantly (refer to Figure 9) . The effect of using 
faster observer poles is shown in Figure 9. In this example 
(taken from [Ref. 9]) the Nyquist plot of the system with full 
state feedback is compared with the plot of the same 
controller with observed state feedback. The optimal filter 
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(a steady state Kalman filter) and an observer with "fast 
dynamics" have the effect of greatly diminishing the stability 
margins of the system. [Ref. 9] 



Figure 9 Nyquist Plots [From Ref. 9] 

Consider the two systems in Figure 10 where both systems 
are minimum phase, servable and controllable. For 
comparison, the controller sections have been highlighted by 
placing them inside the dashed lines. Many of the possible 
transfer functions have been compared and investigated with 
regard to robustness from the connections marked X and XX. It 
has been shown that the closed loop transfer functions from 
r to x are identical along with the loop transfer functions 
(with the loops broken at XX) from u' to u. Yet the loop 
transfer functions (with the loops broken at X) from u" to u' 
are not the same in the two systems unless the relationship 




K[ I + C(SI - A) -1 K] 1 = B [ C (SI - A) ^B] 1 



(4.3) 
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holds. Recall that K is the observer gain which now has 
certain restrictions placed on it from (4.3). This 
relationship is known as the Doyle-Stein (D-S) condition, is 
independent of the control gain and depends only on the open- 
loop characteristics of the observer. When (4.3) holds, the 
observer does not influence the transfer function from r to 
x [Ref. 10]. Finally, note that the loop transfer functions 
are the same only at point XX, which is inside the designed 
controller, but they are not the same at point X, which is the 
control interface to the plant and where natural parameter 
uncertainties will most probably occur. [Ref. 9] 

Doyle and Stein have shown that, if the observer gains are 
parameterized as a function of q (a scalar variable) , then 
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this particular function K(q) must meet certain requirements 
as q-noo with the main one being that 

K (q) /q => BZ (4.4) 
where Z is any non-singular matrix [Ref. 9]. A Kalman filter 
provides for such a function (one of many possibilities) that 
works for all controllable, observable and minimum phase 
systems. This is similar to assuming that the covariance Q of 
the process noise used in the Kalman filter has the special 
form of 

Q(q) = Q 0 + q 2 BWB' (4.5) 
where W is any positive definite symmetric matrix and q is 
the scalar weighting factor defined above. (The value q=0 is 
associated with the nominal steady state Kalman filter gain 
for the assumed process noise Q 0 .) As q-noo, the response tends 
toward the full state feedback response as shown in Figure 11. 
The extra term in (4.5) is equivalent to adding extra process 
noise to the control input of the plant. Normally q is 
increased until a satisfactory compromise between robustness 
and noise performance is obtained. This method gives an easy 
way to balance the requirement between noise rejection and 
stability margins. We have investigated this method with 
regard to the NPS AUV vehicle model and present the results 
in the next section. [Ref. 9] 
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C. SIMULATIONS AND RESULTS 

Following the background development of a Doyle-Stein 
observer (robust observer) that was presented in Chapter II, 
we first need to generate numerical values for the linear 
model of the AUV. Presently, there are many simulation 
languages available, but we used MATLAB. The complex and large 
C program that encompasses the AUV model was implemented as 
a function in MATLAB. This allowed us to try many forms and 
types of controllers and to simulate them in a much easier 
fashion than a conventional language might have afforded. 
Mainly this was because of the rich libraries and simple 
graphics that MATLAB has and the relative ease of changing and 
modifying the algorithms of the controller. 

The AUV model can be fit to an ARX (AutoRegressive) 
discrete time model of 
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A(q'’)y(kT) = B (q 1 ) u (kT) + e(kT) (4.6) 
where A and B are polynomials in the time delay operator q' 1 , 
T is the sampling interval and e is an error sequence 
(possibly colored) . Using this ARX model and a Recursive Least 
Squares (RLS) or Recursive Instrumental Variable (RIV) 
algorithm, the parameters of the model can be estimated. In 
general the RIV methods are better suited to cases when the 
noise is significantly colored (non-white) . In our case, a 
third-order model has been estimated using both RLS and RIV 
with only small differences noted in the numerical results. 
The input to the RLS algorithm is required to be "persistently 
exciting" so that all modes of the model will be excited. 

Since a fourth order model yields very imilar results, 
the choice went to the simpler third order model, further 
confirming the model choice made in Chapter II. Figure 12 is 
a plot of the input, actual output and estimated output. 

The next step was to generate the state feedback and 
estimator gains. MATLAB functions were used extensively in 
this area. For state feedback any pole placement techniques 
will work such as DLQR or PLACE in MATLAB. We used the steady 
state Kalman gains from DLQE for the estimator gains where the 
process noise was adjusted or modified according to equation 
(4.5). Copies of the MATLAB programs used in this thesis can 
be obtained by contacting Professor Roberto Cristi (fourth 
entry in the initial distribution list, phone 408-646-2223) . 
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Figure 12 RLS Model Input/Output Characteristics 
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Figure 13 shows the performance of the controller for 
several choices of robust observer. As outlined in the 
previous section, the robust design yields a family of 
observers (parameterized by the parameter q) which improve 
robustness as q becomes larger. This effect is shown in Figure 
13, where the depth responses obtained with q=0, 10 and 50 are 
shown. Notice the increased stability and reduction of 
oscillations. This is also shown in the input signal which is 
limited to ± 0.4 radians. 

Regardless of the type of controller used (VSC or robust 
observer) , an efficient implementation in hardware requires 
at least a survey of the available AUV architectures. 
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Figure 13 D-S Observer Performance Curves 
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IV. AUV ARCHITECTURE 



A. INTRODUCTION 

AUV systems and their related hardware have become 
increasingly complex in order to satisfy all the levels of 
vehicle control. These levels encompass varying degrees of 
asbtraction from the highest (Artificial Intelligence) to the 
lowest (actuator control algorithms) . Historically, the 
challenge has been to find some sort of hardware and software 
combination to satisfy all the constraints generated by such 
a sophisticated system. The problem is that a general purpose 
microcomputer solution tends to be slow (non-real time) and 
inefficient or non-optimal. The trend is to "tailor" the 
hardware and software to the problem at hand. This 
"application specific" approach is especially germane with 
robotics and AUV systems since they rely heavily on sensory 
feedback where real-time response is a necessity. Most of the 
systems surveyed were set up in some sort of hierarchical 
fashion corresponding to the various levels of "intelligence" 
or abstractions in their mission plans. Each level has its 
own numerical computational requirement; therefore, hardware 
and software selections for each level must be tailored for 
that specific application. This should allow the system to 
be as efficient as possible and to allow the lowest control 
level to run in real-time. This chapter looks at some of the 
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various hardware methods which can help achieve this real-time 
performance requirement at the control-actuator interface and 
yet still allow some flexibility in design which will 
encompass a wide range of applications. 

A . HARDWARE ALTERNATIVES 

The architecture specifications should match or exceed 
the performance requirements of the algorithm which are 
normally formulated in terms of latency (elapsed time from 
when the input is present until the output is ready) and 
throughput rate (given in MIPS or MFLOPS) . There is a 

tendency to concentrate on the throughput performance while 
the latency (i.e, delay) is of critical importance in control 
algorithms (the system stability requirements may allow only 
so much inherent delay) . Control applications normally require 
positional accuracy, concurrent tasks, repeatability, 
robustness and timing constraints (tasks must all be done 
within the sampling period) for the hardware. It seems 
reasonable that the optimal approach would be some combination 
of the various hardware technologies. Some of these 
technologies are: 

• pipelining, 

• RISC (Reduced-Instruction-Set-Computer) , 

• bit-sliced microprocessors, 

• vector processing, 

• DSP (Digitial Signal Processing) chips, 
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• ASIC (Application-Specific Integrated Circuit) chips, 

• systolic arrays, 

• multiple processors, and 

• PLDs (Programmable Logic Devices) . [Ref. 11] 

Pipelining (paralleling the datapath) tries to improve 

the throughput by shortening the clock cycle, but then the 
latency also increases. The concept will also work well at 
higher levels of control or abstraction (concurrent 
programming) and not just at the computational level such as 
a multiplier/accumulator (MAC) . [Ref. 11] 

RISC is a diverse technology and is still subject to much 
debate and company specific design philosophy. The 

performance comes from using less "chip real estate" to encode 
fewer instructions thus leaving more room to add additional 
components to speed up overall program execution. Several 
manufacturers have working RISC microprocessors and their 
throughput performance is impressive. RISC systems have not 
had many commercial applications because twice as much memory 
(versus a 80X86 system) is required and their speed advantage 
has eroded since the newer 80X86 chips are much faster now. 
Even so, AUV control is an area where a dedicated RISC 
architecture might be exploited in a generic controller 
scenario. 

Bit-sliced microprocessors are very fast since they 
generally use an ECL (Emitter-Coupled-Logic) chip set. They 
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can be tailored to an individual application and fit very well 
into fault tolerant systems. Their power hungry ECL 
construction and multi-chip expense limit their general use 
in a space and energy conservative environment such as an AUV, 
but certain performance constraints requiring ultra fast 
processing would benefit from their use. [Ref. 12] 

Systolic arrays and vector processing basically combine 
several arithmetic or processing units in various geometries 
to improve the data stream flow. The general idea is to 
increase performance by some sort of parallel computation. 
Most of the supercomputers use some form of vector processing 
to achieve enormous throughput capability and systolic arrays 
have been used fairly successfully in image processing 
systems. The multi-dimensional nature of an AUV system along 
with the MIMO control problem (with its discrete time 
dynamical difference equations) requires a matrix formulation 
and then a real-time solution which lends itself to this type 
of processing. 

DSP chips have begun to flourish and their performance is 
truly remarkable. Several of the well-established companies 
that market 32-bit floating point processors are: Motorola 
( DSP96002 ) , Texas Instruments (TMS320C30) and AT&T (DSP32C- 
80) . While the state-of-the-art DSP chips are relatively 
expensive, their performance (33 MFLOPS for the T.I. chip), 
onboard memory and on chip I/O make them virtually a single 
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chip solution for almost any control or filter problem. The 
capability of extensive on-chip programming and the ability 
to move blocks of memory onto the chip greatly enhances the 
actual throughput of a large control algorithm. A possible 
pitfall with these chips is that along with the substantial 
performance increases come commensurate leaps in system 
details (debugging programs, etc.). The most effective way 
to minimize this problem and its learning curve is to purchase 
the complete package (C compiler, PC-board, etc.) from the 
manufacturer. Here again it is probably best to stay with 
well-established companies who are marketing their next 
generation chip which already has most of the required 
software available from the previous generation's chip product 
line. The final choice of DSP chip will be a balance between 
performance, cost and supporting software (which is the most 
important item considering the essentially equivalent 
performance from all of the 32-bit chips) . The time and 
effort involved in getting any DSP chip to actually work in 
a system will entail a large portion of the total resources, 
so any product that can shorten or automate this process will 
be money well spent. 

The term ASIC seems to suggest that it is a combination 
of some design methodology, system requirements and choice of 
algorithm which, taken collectively, purport to solve some 
specific control problem. This would allow many of the system 
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overhead functions usually performed by a general purpose 
microprocessor to be eliminated, and hence an overall 
performance increase would be achieved. ASIC devices 
typically are either full custom or semi-custom with the 
latter the most common due to the greatly reduced design 
costs. Semi-custom designs normally use gate arrays or 
standard cells (gates, registers, ALU, etc.) with the former 
the most popular since it generally has a faster design turn 
around time. A natural application for these chips would be 
in an AUV system where navigation, actuator controllers and 
sensor processing could require vast computational power. 
Generally the problems to be addressed are testability of the 
ASIC design and efficient layout (interconnections) of the 
chip to minimize propagation delays. [Ref. 13] 

One of the more exciting areas to explore is that of 
multiple processors. There are various methods and ideas on 
how to implement these distributed computing techniques. One 
of the more important decisions is the one between message 
passing and shared memory. Message passing can have a long 
latency time and a large overhead due to protocol 
requirements. The allure of multiple processors is that the 
type of processor used can be chosen to best complement the 
particular algorithm at hand. How the inter-communications 
between processors is handled is a crucial design decision 
second only to the processor choice itself. [Ref. 11] 
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There are essentially on^.y a few methods of inter- 
communication between processors: RS232 cable, a local or 
system bus, shared memory (dual port RAM, etc.) or a local 
area network (LAN) line (e.g., Ethernet). The LAN tends to 
have a lot of overhead, is serial in nature and, hence, is 
relatively slow. For distant processors that need infrequent 
communications this could be a viable method, but for tightly- 
coupled processors a versatile bus is the only real choice. 
There are several styles and types of buses from which to 
choose; they range from simple passive backplanes to forward- 
looking high-speed 32-bit buses. Issues such as compatibility 
and performance must be investigated and the inevitable 
compromises made. Several types of buses can be employed in 
one system since AUV control systems tend to be hierarchical 
in topology and each layer of control or abstraction can be 
tied together with a particular bus to fulfill a specific 
need. 

The newest 32-bit buses offer the most performance and 
versatility so far. They can block transfer to RAM, do cache 
coherence, autoconf igure (poll boards attached to the bus and 
then adjust the related interface software) and interact with 
the fault-tolerant system (i.e., logically remove the faulted 
board on-line) . Some of the buses can handle larger bus 
widths of up to 256 bits of data (e.g., FUTURE BUS + : IEEE 
standard P896) . The U.S. Navy has decided to base all 
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mission-critical computers on FUTUREBUS. Given this decision 
and the superb performance and adaptability of this bus, it 
would seem a natural choice for future expansions. Of special 
note is the new type of transceiver used by FUTUREBUS called 
BTL (Backplane Transceiver Logic) . This transceiver reduces 
the bus capacitive load and hence increases the bus bandwidth 
to 400 Mbytes/sec (an order of magnitude better than ECL 
transceivers). [Ref. 14] 

Finally, PLDs of which the EPROMs and PALs are examples 
are a relatively cheap way to implement functions or 
processes. PLDs are expandable, universally compatible and, 
by design, tailor-made to the specific application and 
algorithm. These devices are cheap (compared to DSP devices) 
and several versions are in development that have provisions 
to be re-programmed on the fly (EEPROMs , etc.). 

B. OPERATING EXAMPLES 

While there are as many theories and ideas on how an AUV 
architecture should be built as there are AUV manuf actorers , 
it can be quite helpful to investigate a few designs which 
have been built and are actually operational. This section 
will briefly overview three vehicles: 

• ARCS (Autonomous Remotely Controlled Submersible) built 
by International Submarine Engineering Limited (ISE) in 
Canada . 

• EAVE (Experimental Autonomous Vehicle) EAST built by the 
University of New Hampshire. 
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• FS (Free-Swimmer and formally EAVE WEST) built by Naval 

Ocean Systems Center in San Diego, CA. 

Table 1 is a partial list of several AUV designs and their 
ROV manufacturers. Two comprehensive AUV references that 
detail current vehicles are ROV Review 1990 and Undersea 
Vehicles Directly 1990 . [ ROV Review can be purchased after 
January 1990 by contacting Dean Given, PO Box 368, Spring 
Valley, CA 92077, phone 619-660-0402. Undersea Vehicles 
Directory can be purchased after November 1989 by contacting 
Frank Busby in Arlington, VA, phone 704-892-2888.) [Ref. 15] 

The ARCS is a commercially available AUV, so it uses "off- 
the-shelf” hardware and software to reduce costs. The 
modeling of the control system organization followed that of 
a typical naval submarine and a multi-tasking setup was used 
to schedule all the tasks. The specific hardware included a 
16-bit CPU for the multi-tasking, an Intel Multibus and three 
single board computers (two 8086s with 8087s and one 8088) , 
tied together with a common backplane bus and sharing some 
dual port RAM. Interfaces to external equipment were through 
RS-232 serial links. The requirements for a real-time multi- 
tasking operating system led ISE to choose Intel's RMX86 
system. [Refs 16 and 17] 

The University of New Hampshire has been developing and 
enhancing their AUV project since 1977. They are currently 
pursuing their third generation vehicles which stress 
Artificial Intelligence (A. I.) concepts and multi-vehicle 
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cooperation. EAVE EAST III has much greater memory capability 
and processing power than its predecessors in order to serve 
as a testbed for future concepts and evaluations. Their 
computer system is hierarchical with three distinct levels. 
The lowest level handles all the fast (less than one second) 
systems including reading sensors, pre-processing data and 
controlling actuators. This level consists of three 68000 
processors which communicate over RS-232 cable to each other 
and to an interface with the higher levels. This lower level 
can control the vehicle without any higher level assistance. 

The higher level architecture consists of several 68020 
processors on a VME bus and liberal use of dual port memories 
which allow local CPU use as well as communication by other 
CPUs on the bus. Figure 14 shows this multi-bus setup which 
helps to incorporate I/O and data storage into the overall 
system. Their choice of a real-time operating system was pSOS 
which was also tasked with running a symbolic language (LISP) . 
So far their choice of architecture, hierarchy and software 
has proven to be reliable, capable and extensible. [Ref. 18] 

NOSC has a wealth of experience in underwater vehicles 
and their subsystems. Their FS (Free-Swimmer or previously 
known as EAVE WEST) vehicle has been through several upgrades 
in capability. Although it is probably not as advanced as 
some of their current vehicles (which have restricted 
distribution due to assigned mission areas) , FS should be 
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Table I Current AUVs [From Ref. 15] 
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Figure 14 EAVE EAST Architecture [From Ref. 18] 



fairly representative of their architectural design 
philosophy. The NOSC undersea branch is committed to all 
aspects of AUV-related research from using a transputer array 
to simulate an Artificial Neural Network (ANN) to designing 
and testing a Plan Execution System (PES) . [Ref. 19] 

The FS architecture takes a layered approach as pictured 
in Figure 15. Their PES runs on an 80286 CPU and communicates 
with the navigator CPU (8088) and the lower level bus through 
a RS-232 link. The lower level groups several functions 
(using 8088 CPUs) together along with sensors and actuator 
(effector) I/O. Their system requirements lead them to use 
VERTX for the operating system. 
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Figure 15 NOSC FS Architecture 

This is just a sampling of the current operational AUV 
architectural designs. Table 2 lists several other 

organizations and some of their research. [Ref. 12’. 

D. OBSERVATIONS 

The selection of an AUV architecture is a very critical 
process because the system must be expandable and compatible 
enough to meet ever-increasing complex missions needs. The 
design of a generic AUV architecture would entail several key 
items . 

First, the overall system would be hierarchical in nature 
with several layers of control. This aspect would suggest 
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Table II AUV Research Organizations [After Ref. 12] 



VISITED 



PURPOSE/ IMPORTANCE 



AMETEK-STRAZA/SAN DIEGO 
ARL/UNIV. OF TEXAS/AUSTIN 
AUTONET ICS. LOS ANGELES 
DRAPER LABORATORIES. BOSTON 



HYDROPRODUCTS. SAN DIEGO 

EDO WESTERN. SALT LAKE CITY 
GOULD. BALTIMORE 

HONEYWELL MARINE. SEATTLE 

LOCKHEED ADVANCED MARINE 
SYSTEMS, SUNNYVALE 
MARTIN MARIETTA. BALTIMORE 

MARTIN MARIETTA. DENVER 

NOSC, SAN DIEGO 

NOSC, HAWAII 



PERRY OFFSHORE, RIVIERA BEACH 

UNIVERSITY OF UTAH, SALT LAKE CITY 
WESTI NGIIOUSE , ANNAPOLIS 

DEEP SUBMERGENCE LABORATORY OF 
WOODS HOLE OCEANOGRAPHIC 
INSTITUTE. FALMOUTH 



PROVEN CAPABILITIES, TRENDS. WORK SUITES 
OBSTACLE AVOIDANCE SONAR 

NAVIGATION, AOVs. HYDRODYNAMICS. LOW DRAG BODIES 

AUTONOMOUS VEHICLE PROCESSING AND CONTROL. 
NAVIGATION, RELIABILITY ANALYSIS, LIGHTWEIGHT 
STRUCTURES 

PROVEN CAPABILITIES. TRENDS. UNTETHERED 

APPLICATIONS 

SONARS, TV CAMERAS 

ENERGY SUPPLIES. FIBER OPTICS. CABLE HANDLING 
SYSTEMS 

NAVIGATION, ACOUSTIC LINKS. VEHICLE CAPABILITIES. 
TRENCHING, MK-50 TORPEDO, BURIED OBJECT DETECTION 
LOW COST MINE NEUTRALIZATION VEHICLE 



LIGHTWEIGHT STRUCTURES. BONDING. AOVs. IR&D 
PROGRAMS 

INTELLIGENCE PROCESSING. ROBOTICS. AUTONOMOUS LAND 
VEHICLE 

UNMANNED VEHICLE OVERVIEW, AOVs, ACOUSTIC LINKS, 
WORK FUNCTIONS 

DEEP VEHICLES, AUTONOMOUS VEHICLES, FIBER OPTIC 
CABLES, CABLE DEPLOYMENT METHODS. MARINE 
BIOSYSTEMS 

COMMERCIAL VEHICLES. MISSION SPECIFIC VEHICLES. 
WORK PACKAGES, MANIPULATORS 
DEXTEROUS HAND 

PROVEN CAPABILITIES, LAUNCH/RECOVERY SYSTEMS, 
SUBMERSIBLE VEHICLE INTERFACE. SONARS 
DEEP SUBMERGENCE VEHICLES. IMAGING SYSTEMS, 
UNDERWATER SEARCH EXPERIENCE 



several buses linking lower level tasks where they, in turn, 
could also be linked by buses to other hardware or 
subfunctions, etc. This hierarchial design also facilitates 
the implementation of Knowledge Based Systems (KBS) . Choosing 
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the optimum buses becomes extremely important to overall 
system integrity and speed of response. 

Next, the specific type of processor required can be 
chosen to match the particular levels requirements and 
associated algorithms. Whether one uses parallel processing, 
ASIC chips, DSP chips, transputers, or some hybrid combination 
would depend on the cost, speed and expandability required at 
that level. In general, a special function chip could be used 
at the lowest level (dedicated, application specific, stand- 
alone, fast, etc.) to evoke maximum performance while a state- 
of-the-art processor (80X86, RISC, transputer array, etc.) 
would handle the higher level computations which would 
necessarily be constantly changing with the environment and 
mission objectives. 

Finally, the operating system must be chosen to fully 
utilize the hardware's performance and integrate the various 
levels of software. Real-time response, while supporting a 
multi-tasking environment, must be achieved. For future 
growth a symbolic language must also be supported. The 
ability to support limited missions has already been 
demonstrated, but the system level integration and cooperation 
needed to support a high level of Artificial Intelligence (an 
expert system) in an AUV is still being researched. 
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VI . CONCLUSIONS 



A. RESULTS 

Variable structure control is a highly viable control 
methodology for an AUV. Combining an adaptive portion of the 
controller with the VSC method yielded satisfactory 
performance considering the nonlinear nature of the problem. 

Robust observers also have merit in systems with 
uncertainties. In our application, this technique has worked 
reasonably well for a range of speed and model parameter 
changes . 

The ideal AUV architecture would most likely consist of: 

• DSP or ASIC chips at the lowest level tied together with 
a local bus, 

• special purpose processors implementing the higher levels 
and tied together with a smart bus (such as FUTURE BUS + ) , 
and 

• a special purpose optimized real-time operating system to 
implement a multi-tasked hierarchical knowledge-based 
event-driven control system. 

B. FUTURE RESEARCH 

Several areas of AUV control are open to investigation 
and optimization. The use of DSP chips to pre-process inputs 
(sensor fusion) is an important area. The possibility of 
using Neural Networks (and associated parallel processing) to 
handle signal processing and control problems is a very 
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current and far-reaching issue. The appropriate communication 
links (bus type, protocol, network topology, etc.) between 
processors and the various control levels is a critical design 
choice affecting the total system performance. Eventually, 
a fault-tolerant system such as a Structurally Adaptive System 
(SAS) must be looked into to satisfy the reliability issue. 
A Real-Time Knowledge Based System (RKBS) should be researched 
as to the implementation of artificial intelligence into an 
AUV control system. Finally, a MIMO (Multiple-Input-Multiple- 
Output) VSC controller needs to be designed that takes into 
account all of the inputs (depth, speed, etc.), the outputs 
(actuators) and as accurate a model as possible in order to 
maximize the performance. 
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